Size scaling of static friction 
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Sliding friction across a thin soft lubricant film typically occurs by stick-slip, the lubricant fully 
solidifying at stick, yielding and flowing at slip. The static friction force per unit area preceding 
slip is known from molecular dynamics (MD) simulations to decrease with increasing contact area. 
That makes the large-size fate of stick-slip unclear and unknown; its possible vanishing is important 
as it would herald smooth sliding with a dramatic drop of kinetic friction at large size. Here we 
formulate a scaling law of the static friction force, which for a soft lubricant is predicted to decrease 
as f m + A//j4 7 for increasing contact area A, with 7 > 0. Our main finding is that the value of / m , 
controlling the survival of stick-slip at large size, can be evaluated by simulations of comparably 
small size. MD simulations of soft lubricant sliding are presented, which verify this theory. 

PACS numbers: 81.40.Pq; 46.55.+d; 61.72.Hh 



Boundary friction of sliding crystal surfaces across 
atomically thin solid or nearly solid lubricant layers, of 
considerable conceptual and practical importance, also 
constitutes an open physics problem, because the uncer- 
tain occurrence of stick-slip makes the prediction of the 
overall frictional regime - stick-slip or smooth sliding - 
rather uncertain PTII1- While for hard solid lubricants 
the answer is known, namely stick-slip for commensu- 
rate and crystallographically aligned interfaces or smooth 
sliding for lattice mismatched/misaligned interfaces [|[, it 
is not so for soft solid lubricants. The latter, with shear- 
ing occuring inside the lubricant rather than at the sur- 
face/lubricant interface, represent the commonest case, 
realized at room temperature in e.g. commercial machine 
oils confined in between metallic surfaces. The possibility 
of smooth sliding would be especially relevant, because 
of the accompanying large drop of kinetic friction, often 
a very desirable outcome. The crucial controlling quan- 
tity is the magnitude of static friction f s - the maximum 
pulling force reached before slip. So long as f s is finite 
there will be stick-slip; when f s drops to zero, there can 
only be smooth sliding. Realistic molecular dynamics 
(MD) simulations of lubricants confined between atom- 
ically fiat surfaces generally indicate that stick-slip pre- 
vails for soft lubricants, with consequently high kinetic 
friction. However, while in smooth sliding the kinetic 
friction per unit area is essentially size independent, its 
static counterpart f s may decrease with increasing con- 
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tact area A [l],[H. Despite the increased computer power, 
the simulated system sizes are still far too small to 
establish conclusively whether in the limit of mesoscopi- 
cally large size the static friction will remain finite, and 
stick-slip will survive with large kinetic friction, or if it 
will vanish so that smooth sliding and low dynamic fric- 
tion will eventually prevail. The time-honored approach 
borrowed from equilibrium statistical mechanics to this 
type of question is finite-size scaling [8[ . One can for ex- 
ample double repeatedly the size of the simulation cell 
and compare the change in the results with some ana- 
lytically predicted size dependence from theory. Given a 
good scaling prediction, a few simulated sizes are often 
sufficient to establish the large-size limit with reasonable 
accuracy, and in particular whether static friction will 
drop to zero and stick-slip will disappear, or not. 

In this Letter we solve this question, first by deriving a 
size scaling law for static friction, and then showing that 
it fits realistic MD simulations yielding a well denned 
answer. Our end result is that (i) the predicted drop of 
the size-dependent part of the static friction per atom f s 
is inversely proportional to the linear size of the contact 
(i.e. to A 1 / 2 ), but that (ii) its predicted large-size limit 
is nonzero, so that stick-slip will generally survive in soft 
solid boundary lubrication. 

Scaling theory. To start off the theory, we inspect first 
the dynamics of previous MD simulations of sliding over 
soft solid lubricants @, [t| [l(|. These simulations indi- 
cate, very reasonably, that unlike hard lubricants where 
sliding occurs at the interfaces, plastic motion within the 
soft lubricant nucleates typically at some weak point well 
inside the lubricant film, such as a point defect, a dislo- 
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cation, a local incommensurability, etc. (similarly to the 
"weakest-link hypothesis" of fracture mechanics, e.g., see 
Refs. [HI HI] and references therein). The static fric- 
tion force f a (per substrate atom, i.e. f s — a s L x L y /N s , 
where a s is the shear stress, L x and L y are the sides of 
the rectangular simulation cell, and N s is the number 
of substrate surface atoms) depends on the given initial 
(frozen) configuration. For a given size A — L x x L y of the 
simulation cell, different realizations of the initial config- 
uration will give different f s values f s i < f s2 < .. . , where 
a given value f s $ is realized with probability p\ A ^ . Now 
suppose we double the simulation cell. Slip-motion will 
again start at the weakest point wherever it is, in cither 
half of the doubled cell. Assuming that the new (larger) 
contact does not develop new thresholds, the probabil- 
ity that the doubled cell fails at threshold f s < equals the 
sum of the probability that failure occurs precisely at this 
threshold f s j in both halfs plus the probability that in 
one half the threshold is f S i and in the other half it is 
some larger f s f. 



V 



(2A) 



j>i 



(1) 



The factor 2 accounts for the the two symmetric realiza- 
tions of f s j and f s i in the two halves. 

By iteration of Eq. (fTJ), we can find the probability 

p\ AA ^ for larger and larger cell size AA = 2 n A, with 
n = 0,1,... Given the resulting distribution, one can 
calculate the average static threshold for the AA cell by 
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To illustrate this approach, consider the simple instruc- 
tive example where only two thresholds f s i — f m and 
./« 2 = /m + A/ > f m occur, with probabilities pi andp2- 
For the doubled cell, we have four possible thresholds re- 
alizations: (f si, f s i) with probability pf , (f sl , f s 2) with 
probability pip 2 , (/ S 2,/si) with probability p 2 pi, and 
(/s2)/s2) with probability p\. Accordingly, the dou- 
bled cell fails at the lower threshold with probability 
2pip 2 , and at the upper threshold with probabil- 
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The iteration chain is p 2 



(2AA) 



P ( 2 AA) 



a'% and thus p^ A ^ — 1 — a A . Accordingly, the 
average static friction approaches the minimum threshold 
f m exponentially in A: 



a and p 2 
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a < 1. 



P ( 2 AA) 



with solution 



/s(A) — f m = a A Af = e AlnQ A/ . 



(3) 



When we replace the discrete thresholds f S i with a 
more realistic continuous distribution with probability 

Pc ; (fs), the first, quadratic, contribution in Eq. (fTJ) can 
be neglected, and the iteration equation takes the form 
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FIG. 1: (color online) Iteration of Eq. ((4} starting from a saw- 
tooth initial distribution P^ A) (f) = 2 / 0(f) 9(1 - /), where 
9(x) is the Heaviside step function. Inset (a) shows succes- 
sive iterations of Pc (f) for A = 1, 2, 4, 8, . . . Inset (b) dis- 
plays the scaled distribution A~ 7 Pj AA ^ (A~ 7 /) compared to 
the infinite-size solution QD) . The main graph shows the av- 
erage static friction / excess (above the minimum f m ) as a 
function of contact size. 



Figure [TJa) illustrates an example of the numerical iter- 
ation of this equation. Simulations suggest that for large 
A the distribution Pc AA ^ (f s ) tends to approach some uni- 
versal shape, with little dependency on the small-size dis- 
tribution Pc(fs), once it is rescaled appropriately. A 
similar scaling behavior was proven for the strain distri- 
butions of the (related) fiber- bundle models [l3l - fT7| . In 
these models the conditions for the emergence of a critical 
point, i.e. a finite stress in the large-scale limit, were in- 
vestigated under the assumption of a nonzero single-fiber 
breaking probability for arbitrarily small stress. Here in- 
stead we consider that more generally the minimal con- 
tact ("single fiber") distribution of unpinning forces f s 
can start off at a minimum f m which can be nonzero. 
Iteration of Eq. Q guarantees that for any contact size 
AA the distribution Pc AA \f s ) vanishes below the same 

f m as Pc (fs)'- scaling preserves f m . To address the 
scaling of the distribution above f m , it is convenient to 
introduce / = / s — f m . Let us assume that at large A 
the normalized probability distribution scales as 



P^ A \f) = aP^ A \af) 



(5) 



where a > 1 is a constant. By substituting Eq. ([5]) 
into Eq. (j4|), for the large-size distribution g(f) = 
liniA-j.oo A >p( AA ^ (A 1 f), with 7 = log 2 a, we obtain 
the following equation: 



P^ A \fs) = 2P^(f s )x P^(f's)df' s . (4) 



ag(af) = 2g(f) J g(f')df' ) 



(6) 
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or 

Differentiating both sides with respect to /, we get 

K) <?(/) ~ ag(af) g'(f) + 2cf{f) = . (8) 

The solutions of this equation depend on a single fea- 
ture of the distribution g(f), namely its small-/ behav- 
ior. More precisely, assuming that g(f) — J2k=k c kf k > 
with <; > -1, we have that 

a = V with 7 = (1 + ko)- 1 , (9) 

K/)=c^ex P (-^) (10) 

solve Eq. ©. Figure QTb) demonstrates the approach 
of the scaled distributions to the function g(f) = 
f exp(— / 2 /2) (with a = ^/2) obtained by starting off 
with an initial distribution P { C A) '(/) = 2 / 9(f) 9(1 - /), 
i.e. with fco — lj c i = 2. 

The scaling theory makes the following predictions: 
(i) as scaling preserves / m , it is possible to predict the 
minimum threshold f m from an evaluation of Pc A ^ (/) 
at the smallest contact size; (ii) the iteration defined by 
Eq. Q preserves the leading term in the / power expan- 
sion of P c (/) above f m ; (iii) regardless of the overall 
shape of the small-size threshold distribution, for large 
size the distribution acquires the "universal" shape of 
Eq. (Unj); (iv) its width A/( AA ) scales down as an inverse 
power law of A; (v) this power law is dictated uniquely 
by the leading power law with which the arbitrary-size 
threshold distribution behaves for f s near f m ; (v) as A 
increases, the average friction force f s approaches f m ac- 
cording to the law 

/ s (A)-/ m ~ [/ a (A = l)-/ m ]A-T. (11) 

In the example of Fig. [1] this relation yields a mean excess 
static friction scaling as the inverse square root of size 
/ s - f m oc A-Va. 

Simulation. To validate our prediction with MD simu- 
lations, we use our previously developed model pi. l9l. [To| . 
Each of the two substrates is modeled by two atomic 
layers, one rigid and one deformable. In the minimum 
size simulation of contact area A, these substrates are 
composed by 12 x 11 atoms arranged in a square lat- 
tice. The space between the substrates is filled by three 
incomplete layers of lubricant atoms (to prevent crystal- 
lization of the lubricant, we put approximately 90% of 
the atoms which would complete 3 perfect monolayers). 
All atoms interact according to the Lennard- Jones (LJ) 
potential. The strength of the lubricant-lubricant inter- 
action is Vu = 1/9 (in dimensionless "natural" units, n.u., 
defined for example in Refs. @, Q), while the lubricant- 
substrate interaction is much stronger, V s i = 1/3. The 
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FIG. 2: (color online) (a) The time evolution of the spring 
force during a segment of the simulated stick-slip dynamics of 
the 12 x 11 substrate model driven at speed 0.01 n.u., with an 
applied load of 0.1 n.u. per rigid substrate atom; circles mark 
the stick-to-slip transitions where the individual static-friction 
thresholds f 3 i are extracted, (b) The probability distribution 
of the static-friction thresholds Pc (fs) as estimated by a 
histogram of the f s % values. A dotted line marks the mean 
value f s , and an arrow marks the estimated f m . 



equilibrium distance of the LJ lubricant-lubricant inter- 
atomic potential is ra/a s = 3.95/3 (i.e., the solid lu- 
bricant is incommensurate with the substrate). These 
parameters correspond to a soft lubricant. Once the 
thin soft lubricant film is interposed between the slid- 
ers, sliding takes place inside the lubricant (as opposed 
to hard lubricants where the sliding would take place at 
the interfaces) and the film melts during sliding, realiz- 
ing the melting-freezing mechanism of stick-slip motion. 
The bottom substrate is kept fixed, while the rigid top 
slider layer is pressed with a load of 0.1 n.u. per substrate 
atom (representing a pressure in the order 100 MPa if 
the model represented a noble gas solid lubricant be- 
tween metal surfaces) and driven through a spring of 
elastic constant k = 3 x 10 -4 n.u. per atom with a ve- 
locity v. We carry out simulations at driving velocities 
v = 3 x 10~ 3 -f 3 x 10~ 2 n.u. These velocities are suffi- 
ciently small that the system exhibits stick-slip motion, 
as illustrated in Fig.^a). We carry out runs of duration 
exceeding 10 7 n.u. for the smallest-size system (12 x 11), 
representing A = 1. By extracting the "static friction" 
thresholds f s { marked by circles in Fig.^a) in correspon- 
dence to the peaks in f(t) preceding each slip event, we 
obtain a set {f S i} sampling the probability distribution 

P<- A \fs)- We evaluate this distribution by means of a 
histogram of 1063 thresholds obtained during a long run, 
and shown in Fig. [Hb). Although the detailed behav- 
ior near the minimum threshold f m is naturally affected 
by limited statistics, the data are consistent with a dis- 
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FIG. 3: (color online) The average static friction force per 
substrate atom as a function of the system size, as obtained 
from molecular dynamics simulations in the same conditions 
as those of Fig. [2] with several sizes multiple of the 12 x 11 
substrate model which is represented by the first A = 1 point 
in figure. The solid line shows the scaling law, Eq. (|11|) . Inset: 
log- log scale. The dashed line marks f m , estimated by the 
lowest observed slip threshold in all simulations of all sizes. 



tribution staring off at f m ~ 0.0075, with an approxi- 
mately linear slope (fco = 1), which produces an expo- 
nent 7 = 1/2. From the threshold distribution Pc A \f s ) 
we extract the mean value f s = 0.0397, marked by a 
dotted line in Fig.^b). 

Having thus characterized the small-area sliding be- 
havior, we proceed to increase the area in order to track 
the size-induced changes. The cell is successively in- 
creased to A = 2 = 2x 1,3 = 3x1, 4 = 2x2, 6 = 3x2, 
9 = 3 x 3, 12 = 4 x 3, 16 = 4 x 4, 20 = 5 x 4, and 
25 = 5 x 5. The results are presented in Fig. [3l As ex- 
pected, the average static friction decreases with system 
size. As was noted, it would not be feasible to extract a 
large-size limit in the absence of a scaling law. We find 
that scaling law (ITTj) with 7=1/2 fits the simulation re- 
sults with reasonable accuracy. The static friction tends 
to a finite large-size value f m > 0, and therefore stick-slip 
will survive at macroscopic size. 

Discussion. We just arrived at the conclusion that once 
f m starts off nonzero, f s will converge to f m > and 
static friction will not disappear in the large-size limit. 
One should however not jump to the conclusion that once 
the static friction threshold distribution starts from zero, 
f m = 0, the amplitude of stick-slip jumps of f(t) will drop 
to zero, and stick-slip friction will necessarily disappear 



in the limit of large contact area. There are two reasons 
why this is not generally true. The first reason resides in 
statics, and follows from elasticity of the substrate. The 
size X c of a domain that can be considered as rigid and 
slides as a whole, is determined by the elastic correla- 
tion length [l8|, ■ Therefore, the average static friction 
force should reach a plateau for sizes L ~ A 1 ! 2 > A c [2fj| . 
The second reason follows from kinetics. When the slid- 
ing motion starts off at some weak contact site, it may 
either die off, or spread over the whole interface with 
some speed c. This process takes a finite time r ~ L/c. 
If at a given driving velocity v, r is of order or larger 
then the time between successive slips t ss ~ f s (L)/(kv), 
then the local sliding initiated by this weak contact will 
lose its role and effectiveness for all sizes L > Xd, where 
Xd = (c/v) fs(Xd) /k, hence the static friction will satu- 
rate rather than decrease further to f m as predicted by 
the scaling theory Eq. pip . 

This leads us to ask more generally what is lost in the 
scaling approach? We assumed that the doubled cell has 
the same set of thresholds {f S i} as the original one. Yet 
a larger cell may develop new collective excitations, e.g., 
a dislocation loop of a size > L. However, if the original 
cell A is large enough, we may still safely extrapolate 
its distribution to a mesoscopic size, where the master 
equation approach [2l[ is applicable. 

Another situation where a different behaviour is ex- 
pected is the Aubry incommensurate superlubric state, 
which has zero static friction threshold in the infinite 
system [22[. That state, corresponding to f m = in our 
theory, occurs preferentially for hard lubricants whose in- 
terior does not develop a shear band, and which do not 
melt during sliding. 

Experimentally, the scaling behavior predicted here 
could be probed by comparing friction-force microscopy 
realizations with tips of different curvature radii [23l . l24j 
sliding on a surface covered by a lubricant close to its 
melting points, e.g. an octamethylcyclotetrasiloxan or an 
ionic liquid at room temperature, or a noble-gas layer 
at a cryogenic temperature. The technological bottom- 
line lesson, finally, is that stick-slip could be attenuated 
by reducing the smallest threshold force, for example by 
promoting extra defects in the lubricant film by additives 
or other means. 
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